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Observation of long-lived polariton states in semiconductor microcavities across the 

parametric threshold 
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The excitation spectrum around the pump-only stationary state of a polariton optical parametric 
oscillator (OPO) in semiconductor microcavities is investigated by time-resolved photoluminescence. 
The response to a weak pulsed perturbation in the vicinity of the idler mode is directly related to the 
lifetime of the elementary excitations. A dramatic increase of the lifetime is observed for a pump 
intensity approaching and exceeding the OPO threshold. The observations can be explained in 
terms of a critical slowing down of the dynamics upon approaching the threshold and the following 
onset of the soft Goldstone mode. 

PACS numbers: 



One of the most striking consequences of the quantum 
nature of matter is Bose-Einstein condensation, a phase 
transition that does not depend on interactions between 
the microscopic constituents but is driven by their indis- 
tinguishability. A quite general property of Bose-Einstein 
condensed quantum fluids is superfluidity, i.e. the ability 
to flow through a container almost without any friction. 
Since the first discovery in 4-Helium, this effect was in- 
vestigated in a variety of quantum fluids, ranging from 3- 
Helium to ultracold atomic gases (author?) [1]. Inspired 
by recent developments in the study ofphase transitions 
in non equilibrium systems (author?) [2], researchers are 
presently working on the extension of the quantum fluid 
concept to many particle systems whose state is no longer 
determined by a thermodynamical equilibrium condition, 
but rather by a balance between the external driving and 
dissipation. 

In this perspective, exciton-polaritons in semiconduc- 
tor microstructures appear as very promising objects, as 
their mixed excitonic and photonic nature allows for a 
number of remarkable properties, e.g. a very light mass, 
significant interactions, and the possibility of all-optical 
manipulation and diagnostics (author?) [3]. On the 
other hand, the finite lifetime of polaritons requires some 
pumping mechanism to continuously replenish the sys- 
tem, which forbids the establishment of a true thermal 
equilibrium in the gas. Several groups have reported 
the observation of spontaneous coherence in spatially ex- 
tended systems of polaritons far from thermodynamical 
equilibrium with mechanisms that can be interpreted as 
non-equilibrium analogs of Bose-Einstein condensation 
(BEC) (author?) Q. 

In spite of the different pumping schemes (resonant for 
the OPO case (author?) [5|, |y, [O, [8] , non-resonant for the 
so-called polariton BEC case (author?) [9]), a U(l) sym- 
metry is spontaneously broken in all cases and coherence 
is not simply inherited from the pump beam. Away from 
the bistability regime (author?) [10], the OPO transi- 



tion is smooth and shows a critical behavior that closely 
resembles the one of a second-order like phase transi- 
tion (author?) [11]. A typical feature of second order 
phase transitions is the critical slowing down of the decay 
rate of the elementary excitations upon approaching the 
transition from below. When a continuous symmetry is 
spontaneously broken above the transition point, the life- 
time of elementary excitations remains divergent in the 
long wavelength limit in agreement with the Goldstone 
theorem of statistical mechanics (author?) [12Lll3j|. 

In the present Letter, we experimentally investigate 
and theoretically model the lifetime of the elementary 
excitations of a quantum fluid of polaritons as the pump 
intensity is spanned across the threshold in the optical 
parametric oscillator (OPO) configuration. The steady 
state of the system is probed by injecting extra polari- 
tons by means of a weak pulsed beam, and the decay 
time of the response is measured as a function of the 
pump intensity. A dramatic slowing down of the dy- 
namics is observed as the threshold is approached from 
below: close to the threshold, the decay time can be- 
come orders of magnitude longer than the typical life time 
of polaritons, and it remains very long even well above 
the threshold. Good agreement between the experimen- 
tal observations and the theoretical model based on the 
generalized polariton Gross-Pitaevskii equation is found, 
which supports the present theoretical understanding of 
the strongly modified dispersion of the elementary exci- 
tations in presence of a strong pump beam. In particular, 
this observation suggests the possibility of investigating 
the polariton dynamics beyond the limits imposed by the 
intrinsic polariton life time. 

The microcavity sample used in these experiments is a 
GaAs/AlAs-based A/2 cavity with a top (bottom) Bragg 
mirror of 15 (25) Alo.1Gao.9As/AlAs pairs, grown on a 
GaAs substrate. A 20 nm wide GaAs quantum well is 
embedded at the antinode position of the cavity mode. 
The sample is kept at a constant temperature of 10 K. 



The quantum well excitons are in strong coupling with 
the cavity mode, with a Rabi splitting of 2Qr = 4.4 meV. 
All the experiments reported here are performed in the 
resonance region, where the k = cavity-mode has the 
same energy as the exciton state. 

In the OPO configuration, polaritons are coherently 
injected into the microcavity by a pump beam, which 
resonantly populates a polariton mode with a defined 
momentum and energy (pump state). Our pump beam 
is a continuous- wave laser source (TKAI2O3), which ex- 
cites the sample with an incident angle of 10° and has 
a 45 fim spot diameter. Its frequency is chosen close 
to resonance with the lower polariton branch (LPB) in 
a way to inject polaritons with a given wave vector k p 
and energy TojO p . Polariton-polariton binary collisions are 
responsible for the parametric scattering of pump po- 
laritons into a pair of distinct signal and idler modes. 
The efficiency of the parametric processes was optimized 
by tuning the pump at a frequency hwp = 1.5273 eV, 
slightly above the linear-regime LP dispersion (author?) 
[3] (author?) [loj], e L PB(k P ), hw p - e L PB(kp) ~ j(k p ). 
Here 7(fc p ) ~ 0.4 meV is the LPB linewidth at k p , while 
the laser linewidth is around 0.1 meV. An additional 2 
ps-long probe pulse, coming from a different laser source 
with a repetition rate of 82 MHz, is incident on the sam- 
ple with a tunable angle and is focused within the pump 
spot with a smaller, 25 /im, spot diameter. Photolumi- 
nescence (PL) is collected and analyzed by a spectro- 
graph coupled either to a streak- or a conventional C CD- 
camera. The emission in the far field is visualized by 
means of a lens on the Fourier plane (k x ;k y ): a direction 
k y of this plane is selected and energy-resolved by the 
spectrograph, which allows the direct observation on the 
CCD camera of any two-dimensional section (k x ;E) of 
the polariton dispersion. To obtain time-resolved PL im- 
ages a streak camera is used in either a (E; t) or a (k x ; t) 
configuration, which allows for a complete study of the 
dynamics at different energies and momenta. 

Firstly, we have investigated the stationary-state po- 
lariton emission in the absence of the probe. Typical 
energy-momentum emission patterns are shown in FigO] 
for two different values of the CW pump power I p . Al- 
though the pump wavevector lies well outside the /c-space 
region imaged in FigO] a small polariton occupation of 
the LPB bottom still appears as a consequence of incoher- 
ent relaxation processes even at low pump powers (lower 
panel). As the polariton density is very low, the photolu- 
minescence spectrum is concentrated on the linear-regime 
LPB branch. 

At higher pump intensities, polariton-polariton inter- 
actions are able to significantly modify the emission pat- 
tern and, in particular, are responsible for parametric 
processes, where two pump polaritons at k p are trans- 
formed into a pair of signal and idler polaritons of 
wavevectors k s ^, respectively. Thanks to the large size 
of the pump spot, the wavevector is approximately con- 
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Figure 1: Direct experimental observation on the CCD cam- 
era of a two-dimensional section (k x ;E), with k y = 0, of 
the lower polariton dispersion for pump powers I p (a) just 
above (13 mW) and (b) below (12 mW) the OPO threshold 
Ip =12.5 mW. Energy spectra at k x = (dashed line in the 
figure) are depicted on the right, with a magnification of a 
factor 3 for the lower panel of the figure. The PL emission 
is normalised to 1 and plotted in a linear color scale, (color 
online) 



served, k s ~ 2 k p — ki and the unique dispersion of po- 
laritons allows for this process to take place in a res- 
onant way, 2e L p B (k p ) ~ e LPB (k s ) + e L p B (ki). The 
onset of parametric oscillation is clearly visible in the 
emission pattern for pump intensities above the thresh- 
old Ip h = 12.5 mW (upper panel of FigH]): the occu- 
pation of the signal at k s ~ becomes in this case 
very large and the linewidth of the emission in energy 
is substantially reduced as compared to the bare LPB 
linewidth shown in the lower panel. The spectral narrow- 
ing in energy is accompanied by a significant broadening 
of the /c-space emission. The flat shape of the coher- 
ent OPO emission in the (k x ; E) plane is however more 
likely to be a consequence of the peculiar shape of finite- 
size non-equilibrium condensates discussed in (author?) 
[15] rather than an evidence of the diffusive nature of the 
Goldstone mode (author?) [16]. 

To clearly identify the threshold, we have studied the 
energy of the signal emission as a function of pump power 
(full triangles in Figj2]). A smooth and almost linear blue- 
shift of the signal energy is visible at low pump powers, 
while a sudden jump appears for I p just above 12 mW 
due to the onset of parametric oscillation. Such a discon- 
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Figure 2: Energy (full triangles) and decay time (open cir- 
cles) of the signal versus pump power. The red dashed line 
indicates the pump threshold. Inset: time evolution of the 
energy-integrated signal emission A/ s for three pump powers 
(1 mW, 7.5 mW, 12 mW). (color online) 



tinuous behaviour around the threshold was predicted in 
Ref (author?) [Ii|(author?) [10]: below the threshold, 
the energy of the (incoherent) parametric luminescence is 
indeed fixed by the slightly blue-shifted LPB dispersion, 
while above threshold it is determined by a more com- 
plex OPO dynamics that also involves the idler energy. 
Far above threshold, the blue-shift saturates. 

The main object of the present study is the response of 
the system in its stationary state to an additional weak 
(< 0.2 mW) probe pulse that impinges on the sample 
at a large angle of around 20°, i.e. in the vicinity of 
the idler. The evolution of the system in response to 
the probe pulse is monitored by investigating the time- 
and momentum-resolved signal emission and, in partic- 
ular, its decay time. As an example, we have traced 
in the inset of Figj2] the time-evolution of the difference 
A/ s = I s (pump + probe) — I s (pump) between the signal 
emission intensity in the presence and in the absence of 
the probe, respectively, for three different pump inten- 
sity values. To rule out non-linear effects and ensure we 
are in a linear-response regime with respect to the probe 
intensity, we have checked that the physics of interest 
is independent of the probe pulse intensity: while some 
changes remain visible in the short-time dynamics, the 
long-time dynamics simply shows a global rescaling of 
the observed intensity. 

Right after the arrival of the probe pulse, parametric 
scattering of pump polaritons into the k s state is stim- 
ulated by the small population of the new idler polari- 
tons injected by the probe: in the plotted curves, this 
corresponds to a fast switch-on of the A/ s at the probe 
arrival time. The fast decay on a 30 ps scale is then 
followed by a much slower exponential decay on a time 
scale in the 100 ps range, i.e. orders of magnitude longer 
than both the empty cavity decay time (2 ps) and the 
polariton-polariton scattering time (author?) [6||, but 



Figure 3: k x dependence of the signal decay time for two 
different pump powers, 11.5 mW (open circles) and 12.5 mW 
(full circles). The high values of the decay time (> 500 ps) far 
from k x = are due to the integration over all the spot area 
and over all the emitted energies: the intensity profile (not 
shown) for both 11.5 mW and 12.5 mW is strongly peaked 
in the parametric scattering region around the signal state at 
k x — 0, while a much weaker and slower incoherent emission 
coming from higher energy states (e.g. residual excitons) and 
from the border of the spot is responsible for the tails at 
\k x \ > 0.5 /im -1 . (color online) 



still significantly shorter than the probe repetition time 
of approximately 12 ns). 

As easily seen by comparing the three curves shown in 
the inset, the response of the system strongly depends on 
the intensity I p of the cw pump. While the decay time 
of the transient decreases for increasing pump powers 
and eventually goes below the streak-camera resolution of 
^30 ps, the long decay-time significantly increases with 
I p . This latter dependence is summarized by the open 
circles in the main panel of FigEl The decay time shows 
a divergent behavior (author?) [17] for pump powers 
approaching I p = 12.5 mW; for higher pump powers, it 
exceeds the time window of our setup. The parametric 
nature of the enhanced life time is confirmed by the coin- 
cidence of the divergence with the signal emission energy 
jump and the frequency narrowing of the luminescence 
(see Fig. CQ). 

The role of the parametric processes is further evi- 
denced by the momentum-resolved data shown in Fig|3j 
in which the decay time of the different k s components 
of the signal emission is plotted as a function of the 
wavevector k x . The considered wavevector range is cen- 
tered around the value k s where signal emission would 
appear if the pump intensity was above threshold. While 
the decay time is a smooth function of k x for I p well be- 
low the parametric oscillation threshold (open circles), a 
marked peak is apparent in the vicinity of k s (k x ~ 0) 
for pump intensities around and above the threshold (full 
circles). 

A convenient way to interpret the observed slowing- 
down of the response to the probe is to use the co- 



herent polariton model based on a pair of coupled 
Gross-Pitaevskii-like nonlinear wave equations for respec- 
tively the photon and exciton fields i^c x(r, t) (author?) 

EEl, 



. dtpc 






• 7x\ , 



■n R ^j c + g\^x\ 2 H2) 



We follow the dynamics of the system starting from the 
ipx,c(r,t) = vacuum state, c^c(k) is the photon dis- 
persion, while the exciton dispersion is assumed to be 
flat at cox- 7c,x are the decay rates of the cavity-photon 
and the exciton, respectively. The exciton-exciton inter- 
actions are characterized by the nonlinear coupling co- 
efficient g and Qr is the exciton-photon Rabi coupling. 
The driving F(r, t) is proportional to the incident elec- 
tromagnetic field and has to include both the CW pump 
and the pulse probe: once the system has attained its 
stationary-state under the CW-pump only, an additional 
short probe pulse is applied close to resonance with the 
idler. The following response of the system is monitored 
on the most relevant observables, in particular the polari- 
ton distribution in /c-space. For the sake of simplicity, we 
have limited ourselves to the case of a plane- wave pump 
with a well-defined wavevector k p and periodic boundary 
conditions, while the finite spatial size of the probe beam 
is fully taken into account. 

As discussed in Ref. (author?) [10,1128, the approach 
to the OPO threshold from below is signalled by the de- 
cay rate of some mode tending to zero. As a function 
of k, the decay time results then strongly peaked around 
the point where the decay rate is the smallest: the closer 
to threshold, the higher the peak value. These general 
claims are perfectly visible in the numerical result plot- 
ted in FiglH For pump intensities just below I* h , the 
evolution of the ntegrated signal emission AI s (t) after 
the arrival of the probe pulse is characterized by a short 
transient followed by a much slower exponential decay, 
with a time constant that dramatically increases as the 
threshold is approached [panel (a)]. By comparing the 
overall decay time of the integrated AI s (t) [panel (b)] 
with the ^-dependent decay time [panel (c)], it is imme- 
diate to see that the former is determined by the decay 
time of the longest-lived mode, a quantity that increases 
in magnitude and becomes progressively more peaked as 
the threshold is approached. 

This theoretical result is in good agreement with the 
experimental observations for pump intensities in the 
vicinity of the threshold, but some specific attention has 
to be paid at the experimental data for very low pump 
power. In this regime, the theoretical calculations predict 
that the decay time should go back to the bare polari- 
ton lifetime, while a quite long decay time is observed 
in the experiment even well below the oscillation thresh- 
old and at wave vectors far from the signal emission. To 
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Figure 4: Results of numerical calculations. Upper (a) panel: 
time dependence of /c-integrated signal emission A/ s for dif- 
ferent values of I v \ integration is performed over the /c-space 
region surrounding k s . Central (b) panel: decay time of A/ s 
as a function of CW pump intensity. Lower (c) panel: /in- 
dependence of decay rate for different values of I P /ll h = 0.57, 
0.95, 0.99 below threshold. Inset (d): /c-dependence of the 
decay time for I P /Ip h — 1.15, above threshold; the vertical 
dotted line indicates the coherent signal emission wavevector 
k s , at which the decay time is infinite by definition of spon- 
taneous symmetry breaking, (color online) 



explain this behaviour, one can mention the spatial in- 
homogeneity of the system that smoothens the /c-space 
distribution, as well as the presence of residual excitons 
that have accumulated into long-lived states and that re- 
lax down on a long time scale. Clearly, these incoherent 
scattering processes are most important for low pump 
powers, while coherent parametric processes take it over 
as the threshold is approached. 

As we already mentioned, the decay time above thresh- 
old is too long to be quantitative measured with the 
present setup. Numerical simulations do not suffer from 
such a difficulty, and we summarize here the main fea- 
tures that one expects for this regime. As a consequence 
of the spontaneously broken U(l) symmetry, the spec- 
trum of the elementary excitations is characterized by 
the presence of a soft Goldstone mode: as the wavevector 
q = k — k s of the excitation tends to zero, both its fre- 



quency and decay rate tend to zero (author?) [16J]. This 
prediction is confirmed by the numerical results for the 
^-dependent decay rate that we show in FigH^d): once 
again, the peaked structure of the decay rate as a func- 
tion of k — k s is apparent. It should be mentioned that, 
even though the long decay time of the slowest decaying 
mode eventually appears in the very long time behaviour 
of the integrated A/ s (t) , this feature is often masked by 
other, faster decaying branches that contribute in a much 
more substantial way to the observed A/ s signal. 

In conclusion, we have investigated the response of a 
continuously pumped optical parametric oscillator to an 
additional weak pulse in the vicinity of the idler. The 
emission from the signal is resolved in time and its mo- 
mentum distribution is analyzed as a function of the 
pump intensity: when the pump intensity is close or 
above the parametric oscillation threshold, a lifetime or- 
ders of magnitude longer than the intrinsic polariton life- 
time is observed for polariton modes in the vicinity of the 
signal emission. The experimental observations are ex- 
plained in terms of a critical slowing down of the elemen- 
tary excitation dynamics as the threshold is approached, 
and then of the presence of a soft Goldstone mode above 
the threshold. 
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